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Abstract 


We study the spectral expansion of the semigroup of a general stable process killed on the first exit 
from the positive half-line. Starting with the Wiener-Hopf factorization we obtain the q-resolvent 
density for the killed process, from which we derive the spectral expansion of the semigroup via 
the inverse Laplace transform. The eigenfunctions and co-eigenfunctions are given rather explicitly 
in terms of the double sine function and they give rise to a pair of integral transforms which 
generalize the classical Fourier sine transform. Our results provide the first explicit example of a 
spectral expansion of the semigroup of a non-symmetric Levy process killed on the first exit form 


the positive half-line. 
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1 Introduction 

Consider the Brownian motion process W = {Wt}t>Q started from Wq = a; > 0 and let Tq be the hrst 
time when W hits zero. It is well-known (see [3, Appendix 1.3]) that the distribution of Tq can be 
computed as 



( 1 ) 


and that the transition probability density of the process W killed at Tq has spectral representation 



( 2 ) 


To give a functional-analytic point of view, let us introduce the semigroup of the killed process as the 
family of operators {Pt}t>o dehned by 


PtU{x) = E^[u{Wt)l{To>t}], 


(3) 
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and let us define the Fourier sine transform as 


nn(A) = / ^(a;) sin(Aa;)da;, u G 


(4) 


It is well-known that 11 can be extended to an isometry in and then formula (1) gives us a 

spectral representation of the semigroup in 


Pi = ne-*^ n. 


(5) 


The function sin(Aa;) plays an important role in the above expressions. It is the eigenfunction of the 
one-dimensional Laplace operator A (which is simply the second derivative operator, An = n"), with a 
Dirichlet boundary condition at zero. Note that A is the infinitesimal generator of the scaled Brownian 
motion Xt = \/2Wt. 

Brownian motion is a very special process, because it enjoys many useful properties: it is a diffu¬ 
sion process, it is a Levy process (a process with stationary and independent increments), it is also 
a self-similar process. The list of other Markov processes with an explicit spectral representation of 
the semigroup is quite short. We mention here the well-known case of one-dimensional diffusion pro¬ 
cesses [18], branching processes [19], a family of symmetric Levy processes obtained as a time change 
of Brownian motion [15] and the recent work of Patie and Savov on non-selfadjoint Markov semigroups 
| 20 ]. 

Our goal in this paper is to generalize results (1), (2) and (5) to strictly stable Levy processes (which 
we call simply stable processes) killed on the first exit from (0, oo). This class consists of Levy processes 
which satisfy the scaling (or, self-similarity) property: for any c > 0 the process {cXt\t>Q (started from 
Xo = 0) has the same distribution as {Xcat}t>Q- Such processes exist when the stability parameter a 
belongs to the interval (0,2]. When a = 2 we recover the scaled Brownian motion and for a < 2 we 
obtain a two-parameter family of processes with jumps that we will discuss in detail in the next section. 
The scaling property of stable processes is the main reason why they are so popular among researchers 
and why they appear so frequently in various applications originating in Physics, Chemistry and Biology. 
Stable processes (or Levy flights, as they are also known in these fields) occur in modelling such diverse 
phenomena as fluctuations and transport in plasma, turbulent diffusions, seismic series and earthquakes, 
signal processing and financial time series (see review article [7] for a comprehensive list of applications). 


1.1 Stable processes 

A Levy process is usually defined through the characteristic exponent T(z) := — lnE[exp(izAi)]. The 
characteristic exponent of a stable process is given by 


_ |^|ag7ria(l/2-p)sign(2)^ 


2; e 


Here the parameters (a, p) belong to the following set of admissible parameters 

M := {a e (0,1], p G (0,1)} U {a e (1, 2], 1 —l/a<p<l/a}. 


( 6 ) 


(7) 


The parameter a is the same one that appeared in the scaling property discussed above, and it can be 
shown that p = Po(W > 0), which explains why p is called the positivity parameter. Everywhere in this 
paper we denote p = 1 — p, and, more generally we will use the “hat” notation to refer to any objects 
obtained from the dual process X = —X. 
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Our set of admissible parameters excludes processes with one-sided jumps: subordinators or negative 
subordinators (a G (0,1) and p G {1,0}) and the spectrally-negative and spectrally-positive processes 
(a G (1,2) and ap = 1 or q;(1 — p) = 1, respectively). The hrst two cases are not interesting, as the 
processes have monotone paths and the semigroup of the process on the positive half-line is rather simple, 
and the case of spectrally-positive/negative processes is covered separately in Section 3.1. 

Let us consider some special cases of stable processes. When a = 2 we necessarily have p = 1/2 
and the process X in this case is simply the scaled Brownian motion Xt = \/2Wt. When a = 1 we 
can rewrite the characteristic exponent in the form T(z) = sin(7rp)|z| -|- icos(7rp)2;, so that the process 
X can be written as Xt = sin(7rp)Zt — cos(7rp)f, where Zt is the Cauchy process. A stable process 
with p = 1/2 is symmetric (has the same distribution as X) and can be obtained as a subordinated 
Brownian motion. More precisely, let St be an Q;/2-stable subordinator defined by the Laplace transform 
E[exp(—^S't)] = exp(— and independent of the Brownian motion W, then the process {Xt}t>o has 
the same distribution as {\/2Wst}t>o- 

In the general case (except when a G (1, 2}) a stable process is a pure-jump Levy process characterized 
by the density of the Levy measure given by 

u{x) = 

where we have denoted c = r(l -|- a) sm{7iap)/7i and c = r(l -|- a) sin(7ra(l — p))/vr. The Levy measure 
describes the distribution and the intensity of the jumps of the process, and it is connected to the 
characteristic exponent via the Levy Khintchine formula 

T(^) = — f — 1 — izh{x))u{x)dx, 

Jr 

where h{x) is the cutoff function, which is needed to ensure the convergence of the integral. Using the 
cutoff function h{x) = 0 when a < 1 and h{x) = x when a > 1 one can check that the above integral 
representation for T( 2 ) can be evaluated in closed form as given in (6). This construction also works 
when a = 1 and p = 1/2 (in this case we can take the cutoff function h{x) = sin(a;)). 

A stable process X is also a Markov process, and its dynamics can be described by the infinitesimal 
generator L. This operator is dehned for a suitable set of functions u as follows 

Lu{x) = ^lim ^(E^[u{Xt)] - u{x)). 

As we have mentioned above, when a = 2 the infinitesimal generator is simply the one-dimensional 
Laplace operator, Lu = u". In the symmetric case (that is, when p = 1/2) the inhnitesimal generator is 
the fractional Laplace operator L = —(—A)"/^. In the general case, it is a (non-local) integro-differential 
operator having the following form: 

Lu{x) = / {u{x + y) — u{x) — u {x)h{y))i'{y)dy, if a 7 ^ 1, 

Jr 

Lu{x) = — cos{tt p)u'fx) H— sin(7rp) / {u{x + y) — u{x) — u'{x) sin(p))p“^dp, if a = 1. 

^ Jr 

If we denote by X the Fourier transform operator 

Xf{z) = [ e^^^f{x)dz, 

Jr 

then the infinitesimal generator can be represented in a particularly simple form L = —X~^d/{—z)X (see 
[2] [Proposition 9]). In other words, the infinitesimal generator L is a pseudo-differential operator [27] 
and the characteristic exponent T is the symbol of this operator. 
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1.2 Main Results 


Consider a stable process X started from Xq = x > Q and denote by Tq the first exit time from (0, cxd). 
Let {Pt}t>o be the transition semigronp of the process X killed at time Tq: these operators are dehned 
by (3), bnt with X^ instead of Wt- Then Pt are snb-Markov operators, in particnlar, they are contraction 
operators on The transition probability density pt{x,y) is dehned as the integral kernel of the 

operator Pp. 

e d?/, To > f) = Pt{x, y)dy. (8) 

It is known that the transition density exists when X is a stable process (this follows from Theorem 3 
and Example 4 in [4]). 

From the analytical point of view, killing the process on the hrst exit from the positive half-line 
is eqnivalent to imposing a Dirichlet bonndary condition on the negative half-line. Thns, in order to 
hnd the eigenfnnctions of the inhnitesimal generator we need to solve the eqnation Lf{x) = A/(x) for 
X > 0 with the “bonndary condition” f{x) = 0 for a: < 0. The problem of compnting v(t,x) = Ptu{x) 
(for a snitable u) is eqnivalent to solving the following parabolic partial integro-differential eqnation 
dtv(t,x) = L^vityx) snbject to (i) the “bonndary condition” v{t,x) = 0 for x < 0 and (ii) the initial 
condition n(0,x) = u{x). Finally, the transition probability density pt{x,y) is simply the fnndamental 
solution to this partial integro-differential eqnation. 

In order to present onr results, we need to introduce a certain special function 5*2(^) = S 2 {z] a), called 
the double sine function [11, 12]. This function can be dehned via two functional equations 


S2{z 1 ) 


S2{Z) 

2 sin(7rx/a) ’ 


S 2 {z a) 


S2{Z) 
2sin(7rx) ’ 


(9) 


and the normalizing condition S' 2 ((l -|- Q;)/2) = 1. We collect several equivalent dehnitions and various 
properties of the double sine function in Appendix A, here we only mention the following two facts which 
will be used most frequently in this paper: 


(i) The function S 2 {z) is a real meromorphic function having poles at points z G {ma+n : m,n ^ N}; 

(ii) For every 6, c G M have the following asymptotic result 

\S2{b + ialn{P^y)/{ 27 i))S 2 {b - ialn(e''^?/)/( 27 r))| = ^as y ^ Too, 

\y “/"^+0(1 + o(l)), as y-)■ 0+. 

Moreover, the above asymptotic result holds uniformly in b and c on compact subsets of M. 


Now we introduce two functions G and T, which will play an important role in this paper. For 
(a, p) E A and x > 0 we dehne 


G{x) := / e ^P\S2{1 + a + ap/2+ iahi{z)/{2TT))\^<lz (11) 

Jo 

and 

T(x) := sin(xsin(7rp)-I-7rp(l — Q;p)/2)-I-^^^S' 2 (—ap)G(x). (12) 

dvr 

When it will be needed to stress the dependence on the parameters a and p, we will write T(x; a, p) for 
T(x) and G{x;a,p) for G(x). We dehne F and G in a similar way, exchanging the roles of p and p (so 
that T(x) = F{x;a, p)). 
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Note that the integral in (11) converges for a; = 0 (and thus for all a; > 0): this is easy to establish 
using (10). Another important observation is that the function G is completely monotone, that is 
(—l)"'G^"'^(a;) > 0 for all n G Z’*' and a; > 0. Finally, we note that the function F is bounded on (0, oo) 
if P > 1/2, and it grows exponentially if p < 1/2 (with a similar result for F and p). 

The following theorem is our first main result in this paper: here we generalize formulas (1) and (2) 
which hold in the Brownian motion case. 

Theorem 1. Let X be a stable process defined by parameters {a, p) G A. 

(i) If a > 1 or p >1/2 then for all x > 0 

P^(To >t) = ^S2{ap) f e-‘^“F(Ax)A-MA. (13) 

TT Jo 

(ii) If a > 1 or p = 1/2 then for all t,x,y > 0 we have 

2 

pt{x,y) = - / e-*"“F(Aa:)F(Ap)dA. (14) 

^ Jo 

In the symmetric case Theorem 1 was established in [15] (see Example 6.1), and in the non-symmetric 
case the eigenfunctions F{x) were computed in [16] (in both papers the function G{x) was given in an 
equivalent integral form). 

Our next goal is to study the transition semigroup Ft and to establish an analogue of identity (5). 
In the case of non-symmetric stable processes the situation is bound to be more complicated than in the 
symmetric case, since one of the functions F, F is exponentially increasing. In order to properly dehne 
the operators which diagonalize the transition semigroup, we hrst need to introduce a suitable space of 
test functions. 

Definition 1. Set ( := | min(l, l/a). Let us denote by Ha the set of functions u{x) satisfying the 
following conditions 

(i) u{x) is analytic in the sector | arg(a:)| < ( and it takes real values on (0, oo); 

(ii) For every e G (0,0 there exists 5 = 5{e) > 0 such that \u{x)\ = OdxO^I'^l) as |a;| —)■ oo and 

\u{x) \ = 0(1) as |x| —)■ 0 (uniformly in the sector | arg(x)| < ( — e). 

The sets Ha are non-empty: for example, u{x) = (1 -|- x)~^ G Hi and v{x) = exp(—x") G Hy for 

(3 > a > 1 (but not for /3 < a or a < 1). It is clear that Ha = "Hi for a < 1 and Ha C Hy for 

1 < a < fi < 2. The following properties follow easily from the definition: if u and v belong to Ha then 
the same is true for functions 


(i) u{x)v{x)] 

(ii) au{x) + bv{x) 

(iii) u{ax + b) 

(iv) x“e-''"“M(x) 


for all a, 6 G M; 
for all a > 0,6 > 0; 
for all a > 0, 6 > 0 


We also record here the following important property: 

(v) The restrictions u\^+ of m G Ha are dense in 
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The above property is easy to prove. Indeed, assuming that there exists v G which is orthogonal 

to all u G T-iai we obtain v(x)(l + x)~°‘^dx = 0 for all a > 0. Substituting for (1 + x)® and writing 
w{y) = n(x)^, we see that w{y)e~°''^dy = 0. Therefore, w{y) = 0 for almost all y > 1, and so 
v{x) = 0 for almost all x > 0. 

Next, we define the following integral operators, which generalize Fourier sine transform (4): 

poo poo 

nM(A) = '\/2/7r / F(Ax)M(x)dx, flM(A) = \p2ji / F(Ax)M(x)dx. (15) 

We denote by Pt the transition semigroup of the dual process X, killed on the hrst exit from (0,cxd). 
According to Hunt’s switching identity (see [2, Theorem 5]), the kernel of the dual semigroup is given 
by Pt{x,y) = pt{y,x), thus Pt is simply the adjoint operator of Pt- 
The next theorem is our second main result. 

Theorem 2. Assume that [a, p) G A and p > 1/2. 

(i) A can he extended to a bounded self-adjoint operator H : L^(M’'') i—>• and II : Pa t T^(M’'') 

is a symmetric operator such that IlPa is dense in L^(M+). For all u G Pa we have 

HIIm = u. (16) 


(a) For u G Pa and t > 0 we have 

Ptu = ne-*^“nM, (17) 

and 

UPtflu = (18) 

Note that in the symmetric case p = 1/2 we have H = H, and then formula (16) implies that H is an 
isometry on L^(M+). This result was established in greater generality in [15]. 

Corollary 1. Assume that {a,p) & A, p> 1/2 and A > 0. 

(i) The functions ux{x) := F{Xx) are the eigenfunctions of the semigroup Pt, that is PtU\ = 

(a) The functions F{Xx) are the eigenfunctions of the dual semigroup Pt, in the sense that for all 
V G Pa we have Pttiv = Ile“*^“x. 

Our third main result gives Laplace and Mellin transforms of the eigenfunctions. 

Theorem 3. Assume that {a, p) G A. 

(i) For z > max(0, cos(7rp)) 

e~^^F(^x)dx = ^^S 2 {ap) + a/2 + ap/2 + ialn(^)/(27r))|^. (19) 

The corresponding result for F{x) can be obtained from the above formula by replacing p ^ p. 

(a) If p> 1/2 then for Re{z) G (—a/5,0) 



x^ ^F(x)dx 


r(z)^2(^) 

2S2{ap + z) 


( 20 ) 
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The two formulas (19) and (20) can lead to many useful explicit results about the eigenfunctions. For 
example, when a is irrational, one could use the methods of [9, 10, 14] and obtain complete asymptotic 
expansions of F{x) as a; —)■ 0+ or x —)■ +cxo, as well as convergent power series representations. We leave 
this investigation for future work, and in this paper we only include the discussion of Doney classes in 
Section 3: as we will see, in this case all expressions involving the double-sine function reduce to simple 
hnite products. 

The paper is organized as follows. In Section 2.1 we review some results from the fluctuation theory of 
Levy processes that will be required later on, and we also outline the plan for proving our main results. In 
Section 2.2 we study the Wiener-Hopf factors of stable processes and establish some preliminary results. 
In Sections 2.3, 2.4 and 2.5 we prove Theorems 3, 1 and 2, respectively. In Section 3 we simplify all 
of our formulas in the case when the process belongs to one of Doney classes, while the special case of 
spectrally one-sided processes is treated in Section 3.1. Finally, in Section 4 we present some concluding 
remarks. 


2 Proofs of Theorems 1, 2 and 3 

2.1 The plan for proving our main results 

Our proofs are based on the Wiener-Hopf factorization, which is a key result in the fluctuation theory 
of Levy processes [2, 17]. We present this result here for the sake of completeness. We start with a Levy 
process X (a one-dimensional process with stationary and independent increments) and we denote by 

Xt = sup{W : 0 < s < t}, 2Lt = infjXs : 0 < s < f} 

the running supremum/inhmum processes. We will denote by e(g) an exponential random variable with 
expected value 1/q, and assume that e(g) is independent of the process X. The following result is 
contained in Theorem 6.15(i) and identity (6.28) in [17]: 

Theorem 4 (The Wiener-Hopf factorization). Let X be a Levy process started from zero. Then for 
q > 0 

(i) the random variable We(g) is independent of X^^q) — Xe(g); 

(a) the random variables We(g) — ^e{q) o,nd 2Le{q) have the same distribution. 

Let us explain how we will use the above Wiener-Hopf factorization result in order to obtain in¬ 
formation about the transition probability density of the process killed on the hrst exit from (0,cxd). 
Assume that Xq = 0 and denote by fx{x) and /x(x) the probability density functions of Xe(i) and 
~^e(i)) respectively (we will prove later that these densities exist and express them in terms of the 
double sine function S' 2 ). Due to the scaling property of stable processes, the random variable Xe(g) has 
density with a similar result for X^i^q)- Let us denote by Hq{x,y,z) the joint density of 

(Xe(g), Xg(-^)) for the process started at Xq = x\ 

IPx(Xe(g) e d?/,Xg(g) e dz) = Hq{x,y,z)dydz, where z < mm{x,y). 

According to the Theorem 4, the pair {Xe(^qp 2Le{q)) has the same distribution as {x + Sq — Iq,x — Iq), 
where Sq and Iq are independent random variables having the same distribution as X^^q^ — x and x — X^ ^^y 
respectively. Therefore the function Hq(x,y,z) can be written as follows 

Hq{x, y, z) = q^^^fxiiv - z)q^^'")fx{{x - z)q^^^). (21) 
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Let us denote by ht{x,y,z) the joint density of {Xt,2Lt) f^e process started from Xq = x, that is 

£ dy, G d^;) = ht{x,y, z)dydz, where 0 < 2 ; < min(a;, ?/). 

This function is related to the transition probability density and the semigroup of the killed process via 
the following identities 


Pt{x,y) = 


min(a:,y) 


ht{x,y,z)dz, Ptu{x) = 


ht{x,y,z)u{y)dy 


dz, 


( 22 ) 


Jo Jo IJz 

where we have assumed that u is bounded on (0, cx)). At the same time, by conditioning on the random 
variable e{q) (which has exponential distribution with the density ge“'^*l{i>o}) we see that 


qe ‘^^ht{x,y,z)dt = Hg{x,y,z), (23) 

and therefore the function q 1 —>• q~^Hq{x,y, z) is the Laplace transform of f 1 —)■ ht{x,y,z). Our plan for 
proving Theorem l(ii) and 2(ii) is to invert the Laplace transform in (23) and to obtain 



ht{x, y,z) = ^ j Vx((l/ “ z)q^^'^)fx{{x - z)g^/")e^Mg, (24) 

JiR 

with fx{x) and f^{x) expressed in terms of the double sine function, and then deform the contour of 
integration, so that the vertical line iM is transformed into Hankel’s contour (beginning at — 00 , going 
around 0 in counter-clockwise direction and ending at — 00 , see Lemma 2 below). As we will see, a lot 
of effort is required for justifying this transformation of the contour of integration. 

Remark 1. The function Hq{x,y,z) is closely related to resolvent operators, defined as 

POO 

Rgu{x) = / e~'^^Ptu{x)dt = g“^Ea;[M(Xe(g))l{x^ >0}]. 

Jo 

It is clear from (22) that Rg is an integral operator with the kernel 

prmn{x,y) 

rq{x,y)=q~^ Hq{x,y,z)dz. 

Jo 


2.2 Some preliminary results 

In the next theorem we identify the Wiener-Hopf factors for a stable process X, which are defined as the 
Laplace transform of positive random variables Xe(i) and — 


(j){z) = (t){z) = E[e^—«(!)], Re(^) > 0. (25) 

Theorem 5. Assume that [a, p) G A. For Re( 2 ;) > 0 

0 ( 2 ;) =2;““^^^S'2(1/2 -I- a/2 + ap/2 + la ln(2;)/(27r)) (26) 

X S'2(1/2 + a/2 + ap/2 — iaTn(2:)/(27r)), 


and (f){z) can be obtained from the above formula by replacing p p- 




The above result was established in [14] in the case when a 7 ^ 1 (see Theorem 4 in [14] and formula 
(A.2) below). Here we present a much simpler proof which also covers the case a = 1. 

Proof of Theorem 5: The proof is based on the Wiener-Hopf factorization result, Theorem 4. Writing 
Xe(i) = dfe(i) + (-Ae(i) “ -^e(i)) and usiug properties (i) and (ii) we obtain the following factorization 


1 

1 + 


E[ei^^e{i)] = E[e'^^-=(1)] X E[e'^^«(i)] = 0(-H) x 0(H), z e M. 


(27) 


This is a classical Riemann-Hilbert problem: we need to hnd a function 0(—iz) analytic in the upper 
half-plane Im( 2 ;) > 0 and continuous in the closed upper half-plane and another function 0 (i 2 ;) having 
the same properties but in the lower half-plane Im( 2 ) < 0 which meet at the boundary Im(x) = 0 as 
prescribed by (27). 

Let us dehne by f{z) the function in the right-hand side of (26), and by f{z) the same function, but 
with p replaced by p. First we will verify that the functions / and / satisfy (27). Assume that z > 0 
and let us denote w = aln(z)/(27ri). Then 

f{-iz)f{iz) = -f- a/4 + ap/2 + w)S2{l/2 + 3a/4 + ap/2 - w) 

^ ^-«(i-p)/2g-™(i-p)/4^^^^/2 + 5a/4 - ap/2 + w)S2{l/2 + 3a/4 - ap/2 - w) 

= + 3a/4 + ap/2 - w)S2{l/2 + 5a/4 - ap/2 + w). 


where we have used identity (A. 7) in the form 


52 ( 1/2 -|- a/4 -|- ap/2 + w)S2{l/2 -|- 3a/4 — ap/2 — w) = 1. 
Next, according to the second functional equation in (9) we have 

^ 2 ( 1/2 + a/4 - ap/2 + w) 


^ 2 ( 1/2 -|- 5a/4 — ap/2 -|- tc) = 


2 sin(7r(l/2 -|- a/4 — ap/2 + w)) 
Using the above result and identity (A.7) in the form 

^ 2 ( 1/2 -|- a/4 — ap/2 + w)S2{l/2 + 3a/4 -|- ap/2 — w) = 1 


we obtain 


1 _ 1 

^o/2g7ria(l/2-p)/2 _|_ ^-a/2g7rio(p-l/2)/2 ]_ _|_ ^ag7ria(l/2-p) ]_ _j_ V[/(^) ’ 

Thus we have verihed that that the functions / and / satisfy (27) for z > 0, and the proof for x < 0 
follows by taking the complex conjugate in the above equation. 

To prove that our candidate solutions / and / are in fact the Wiener-Hopf factors, we need to 
apply a certain uniqueness argument. Using the fact that 52 (^) is analytic and zero-free in the strip 
0 < Re (a) < 1 -|- a we check that the functions / and / are analytic and zero-free in the half-plane 
Re(z) > 0. The asymptotic result (10) easily gives us z~^ ln(/(z)) —)■ 0 as z —>■ 00 uniformly in the half¬ 
plane Re(z) > 0 (with a similar result for /). To hnish the proof we only need to apply the uniqueness 
result given in [13, Theorem 1(f)] and conclude that 0 = /. □ 


f{-iz)f{iz) = 


,-a/2„7ria(p-l/2)/2 


2sin(7r(l/2 -|- a/4 — ap/2 + w)) 

^-a/2g7rio(p-l/2)/2 
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Proposition 1. Assume that (a,p) G A. The random variable -^e(i) 'Is a mixture of exponentials and 
its density is given by 


where 


fxi^) — I ^ X > 0 


pi{u) := — sin(7rap)M"^/^|S'2(l/2 + a + ap /2 + ia ln(M)/(27r))p. 

TT 


(28) 


(29) 


The function fx{.x) (the density of can be obtained from above equations by replacing p p- 

Proof. According to [23, Theorem 2], the random variable ^e(i) is a mixture of exponentials, therefore 
we can write its density in the form (28) with some positive measure p{du). Applying Fubini’s Theorem 
we see that 

^ p{du) 


(j){z) = 


e ^''fxi^)dx = 


u + z 


thus the Wiener-Hopf factor 0(z) given by (26) is the Stieltjes transform of the measure p{du). Using 
the well-known result on the inversion of Stieltjes transform we conclude that the measure p{du) has a 
density p{u) which can be found via 


p{u) = - \m[(j){e^A)]. 

TT 

From formula (26) we hud (as before, denoting w = ia\n{u)/{2x)) 

(j){e^'u) = M“"^/^e“™^'^^S' 2 (l /2 + a + ap /2 + w)S 2 {l /2 + ap /2 — w). 
Applying the second functional equation in (9) we check that 


(30) 


S'2(1/2 -|- ap/2 — w) = 82 ( 1/2 -1- a -|- ap/2 — t(;)2 sin(7r(l/2 -|- ap/2 — w)) 

= 82 ( 1/2 + a + ap /2 - w) -h . 

Combining the above two formulas we hnally obtain 

(j)(e^A) = 1 ^ 2(172 + a + ap /2 + ialn(M)/(27r))|2 ^ 

and now the desired result (29) follows easily from (30) and the above equation. □ 


Let us denote by TZ the Riemann surface of the logarithm function. In what follows we will often 
consider functions defined on TZ or on sectors in TZ. In particular, we have ln(e''^z) = ic -|- ln(z) for all 
c G M and z E TZ. Next, we state two simple lemmas, which will be used often in this paper. Both of 
these results are well-known and can be easily established by Cauchy Residue Theorem. 

Lemma 1 (Rotating the contour of integration). Assume that the function f(z), z eTZ, is analytic in 
the sector —e < a.Tg(z) < b + e for some b > 0 and e > 0, except for a finite number of poles at points 
z = Zi lying in the sector 0 < arg( 2 ;) < b. Assume also that for some 6 > 0 we have f(z) = 0(\z\~^~^^) 
as \z\ — )■ O’*" and f(z) = 0(\z\~^~^) as \z\ — )■ -|-cxd, uniformly in the sector 0 < aig(z) < b. Then 

POO POO 

/ f(z)dz = e* / f(e^’’z)dz + 27ri Res(f(z) : z = Zi). 

Jo Jo 
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Lemma 2 (Hankel’s contour of integration). Assume that the function f{z), z &7l satisfies f{z) = f{z) 
and is analytic in the domain 7r/2 — e < | arg(^)| < tt + e for some e > 0, except for a finite number 
of poles at z = Zi in the sector 7r/2 < arg(^) < vr and the corresponding poles z = Zi in the sector 
—71 < a.Tg{z) < —7r/2. Assume also that for some 6 > 0 we have f{z) = as \z\ —)■ O'*" and 

f{z) = 0{\z\~^~^) as l^l — )■ + 00 , uniformly in the sector ti/2 < a,Tg{z) < tt. Then 

1 r 1 r°° r . 1 

— / f{z)dz = -/ Im[/(e’"‘^)]d 2 ; + 2Re y'Res(/(^) : z = Zi) . 

27ri i* 71 Jo L “ J 

The following elementary result will be useful in deriving various estimates. 

Lemma 3. Let fi{x) := min(x“,a;^) and f 2 {x) := min(a ;'^,where we assume that a > b, c > d, 
a + c + 1 > 0 and 6 + d + 1 < 0. Define f := min(a + 1, —d) and rj := max(6 + 1, — c). Then ^ > rj and 
there exists C = C{a, b, c,d) >0 such that fi{x)f 2 {x/X)dx < C min(A^, A'^) for all A > 0. 

The proof of the above result is very simple - one only needs to evaluate the integral fi(x)f 2 (x/A)dx, 
distinguishing the two cases A > 1 and A < 1. We leave all the details to the reader. 


2.3 Proof of Theorem 3 


Proof of Theorem 3(i): We denote 

/(^) := + a/2 + ap/2 + ialn(z)/(27r)) (31) 

X S'2(l + a/2 + ap/2 — ialn(^)/(27r)). 


We consider f{z) as a function on the Riemann surface JZ. In order to prove part (i) of Theorem 3, we 
will establish an equivalent statement 

FM wy‘dz, (32) 

27ri Jc+iR 

where c > max(0, cos(7rp)). 

Let us denote p = 27rmin(l, 1/a) + pvr (it is easy to check that r/ > vr for all (a,p) G A). Given the 
fact that the double sine function S 2 {z) has poles at points {m + na : m, n > 1} and the pole at 1 + a 
is simple (see Appendix A) we see that the function f{z) is analytic in the sector | arg( 2 :)| < p, except 
for two simple poles at points z± = exp(±7rip). Let us compute the residues at these poles. From the 
hrst functional equation in (9) we hnd that 


S'2(l + a/2 + ap/2 — ia ln(2;)/(27r)) 


S' 2 (a /2 + a()/2 — iQTn(2;)/(27r)) 
2sin(7r(l/2 + p/2 — iln(2:)/(27r)) ’ 


We also check that 

—2sin(7r(l/2 + p/2 — iln(2;)/(27r)) 

Cl >2 

Combining the above two formulas with (31) we obtain 


2 = 2 :+ 


i 


Res(/ : ^ = z+) 


a/q; 

-i-^52(ap) exp(-7rip(l + ap)/2 + 7rip)S'2(a)S'2(l + ap) 
-^exp(7rip(l - ap)/2), 
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\-ap-l\ 


(33) 


where we have also used formulas (A.6) and (A.7). 

Note that 

\f{z)\ < C{a,p) X min(l, |x| 

as \z\ —)■ cx) or |z| —0 uniformly in the sector | arg^l < rj (this upper bound follows from (10) and (31)). 
In particular, the integral in (31) converges, and we can shift the contour of integration in (32) so that 
c + iM H->• iM and then pass to Hankel’s contour of integration (see Lemma 2) and obtain 


— / f{z)e^^dx = -/ e-"^Im[/(xe")]dx + 2Re 


27ri 


' c+iR 


TT 


-- exp(e""^^a; + 7rip(l - ap)/2)) 


_ gCOs(7rp)tc 


sin(a;sin(7rp) + 7rp(l — ap)/2) -/ e ^^lm.[f{ze'^^)]dz. 


TT 


(34) 


Using the definition of f{z) in (31) we check that 


f{ze'"^) = ^^^S 2 {ap)e ^'^^S' 2 (l + a + a;p/2 — ialn(2;)/(27r))S'2(l + q;/)/ 2 + ialn(2:)/(27r)). 

From the second functional equation in (9) it follows that 

52(1 + ap/2 + ialn(x)/( 27 r)) = i [e™/ 5 / 2 ^-a /2 _ g-™p/2^a/2j ^ ^ ^ + ia ln(x)/( 27 r)). 

The above two equations give us 


Im[/(ze’^^)] = sin(7rQ(/))52(ap)2:"^'^^ ^A|5'2(l-)_ q/- f Q/p/2 + iQ(ln(z)/(27r))p. 

Applying the the second functional equation in (9) we check that 


^sin(7rQ;p)52(Q;p) = — ^S2{-ap). 


Combining the above two results with (34) we obtain formula (32). 


□ 


Lemma 4. Assume that (a,p) G A. Then F{x) = 0{x°‘p) as x ^ . 

Proof. Let us consider the function G{x) dehned by (11). Using (10) we check that the integrand in (11) 
satisfies 


zGp-dP\S.,(l + a + ap/2 + ialn(z)/(27r))|2 


z“(l + o(l)), as z —)■ 0+, 

.^-"^-1(1 + o(l)), as z —)■+ 00 . 


(35) 


This result combined with ( 11 ) implies that G( 0 ’'') < +(X) and G"( 0 ^) = —oo. We conclude that the 
function F{x) is monotone in some interval ( 0 ,e). Using ( 10 ) we check that the function f{z) (the 
Laplace transform of F{x) given by ( 19 )) satishes f{z) = Cz~'^^~^{l + o(l)) as z ^ +cxd, for some 
constant C = C{a, p) > 0 . Applying Karamata’s Tauberian theorem followed by the Monotone Density 
Theorem we conclude that F{x) = Cx°‘^{l + o(l)) as a; —)■ O’*". □ 
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Proof of Theorem 3(ii): Let us again denote by f{z) the Laplace transform of F{x), given by (19). 
From formulas (10) and (19) we check that the upper bound (33) is true for all 2 ; G (0, cxd). At the same 
time, when p > 1/2 the function F{x) is bounded and according to Lemma 4 it satisfies F{x) = 0{x‘^^) 
as X —)■ O’*". Therefore we can apply Fubini’s Theorem and conclude that 

poo poo 

/ f{z)z~^dz = — s) / x'^~^F{x)dx, 

Jo Jo 

where both integrals converge absolutely for Re(s) G (—op, 0). The integral identity (A. 11) implies 

/ z-J{z)dz = — - -ap < Re(s) < 1. (36) 

Jo 02(1 — s)o2(a;p + s) 

Formula (20) follows from the above two equations by using (A. 7), the second functional equation in (9) 
and applying the refiection formula for the gamma function. □ 


2.4 Proof of Theorem 1 

In the next lemma we collect some properties of the function p{u) defined by (29). 

Lemma 5. Assume that {a,p) G A. 

(i) The function p{u) is analytic in the sector \ arg(u)| < 7r(l/a + p), u G IZ, except for two simple 
poles at points u± = exp(±7ri(l/a — p)). The residues at these poles are given by 


Res(p(u) U = U±) = ‘^2(«P) gq.^i(app/2+3p/2-l/«) 


2'^^^/a 


(37) 


(a) Denote C-.= a p)/tt. Then 

iM“)i = 


fH“(C' + o(i)), 

as M 

\|u|-“^(C' + o(l)), 

as M 


(38) 


uniformly in the sector \ arg(M)| < 7r(l/a + p). 

Proof. Let us rewrite the expression in (29) in the form 

p(m) = S2{1/2 + a + ap/2 + iQTn(M)/(27r))S'2(l/2 + a + ap/2 — ia!ln(M)/(27r)). (39) 


The double sine function 5*2(^) is a meromorphic function which has poles at points {m + an : m,n > 1}, 
and the pole at z = 1 + a is simple (see Appendix A). This implies that the function p{u) is analytic 
in the sector | arg(M)| < p := 7i/a — vrp + 27rmin(l, 1/a), except for the two simple poles at u = u±. 
It is easy to see that p > 7r/a + up. Let us compute the residue at u = u_. In this case the function 
S' 2 ( 1/2 + a + ap/2 — ia\n.{u)/{/2 ti)) is analytic in the neighborhood of u = and the pole comes from 
the other factor S'2(l/2 + a + ap/2 + iaTn(M)/(27r)), which we transform with the help of (9) as follows 


^ 2 ( 1/2 + a + ap/2 + ia \n.{u)/{/2 ti)) 


S'2(1/2 + ap/2 + ia ln(u)/(27r)) 
2sin(7r(l/2 + ap/2 + ia ln(M)/(27r)) 
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Note that 


thus 


d 

du 


[2 sin(7r(l/2 + ap/2 + ia ln(M)/(27r))] 


U=U — 


a 

m_ ’ 


Resf/ifu) : u = uJ) = —— x -— x (m_)"^/^S' 2(1/2 + a + ap/2 — ialniuJ)/{2 tt)) 

71 a 

X S'2(1/2 + ap/2 + ialn(M_)/(27r)) 

^ sin(7rap) ^,i/2-^i(i/a-p)(i+ap/2)g^^^ ^ ap)52(l) 
vra 

_ 7ri(app/2+3p/2—l/g) 

“ 271^ 


where in the last step we have used the functional equation (9) and formula (A.6). Thus we have proved 
(37) for u = U-, and the result for u = follows by taking the complex conjugate. 

The result in item (ii) follows at once from (10) and (39). □ 


Lemma 6. Assume that (a,p) G A. 

(i) The function fx{x) can be extended to an analytic function in the sector \ arg(a;)| < 7r(l/a + p), 

X e 77. 

(ii) There exists a constant C = C{a,p) such that for all x & IZ in the sector \ arg(x)| < 7r(l/a + p) 
we have 


I „ , . I < f Cmin(|x|"^ Mx| " ^), if | arg(x)| < 7 r(l/a - p), 

^ |Cmin(|x|"^“^, |x|“"“^) + if | arg(x)| > 7 r(l/a — p), 

where we have denoted f := cos(| arg(x)| — 7 r(l/a — p)). 

(in) For x > 0 

= As2(1 + ap) I^F(x) + e"'’F'(a:)) . 


(40) 


(41) 


The corresponding results for /x(x) can be obtained by changing p ^ p and F F. 

Proof. Let us hrst prove part (i). We start with the equation (28), which dehnes fx{x) as an analytic 
function in the right half-plane Re(x) > 0. Assume that x lies in the hrst quadrant, that is arg(x) G 
(0,7r/2). Choose any fd G (0,7r/2) such that fd ^ 7r(l/a — p). Applying Lemmas 1 and 5, we rotate the 
contour of integration in (28) by angle jd in the clockwise direction, so that M+ h-)■ Taking into 

account the pole of p{u) ai u = U- (which will lie in the sector —/d < axg{u) < 0 if /7 > 7r(l/a — p)) we 
obtain the following identity 

poo 

fj^{x) = e~'^ - 27ri X Res(p(M) : n = n_)e"(“-^^l{; 3 >^(i/„_p)}. (42) 

Jo 

Note that the integral in (42) converges if Re(e“‘^x) > 0, thus we have obtained an analytic continuation 
of/(x) into the half-plane arg(x) G {(d — 7 t/ 2, jd + 7 t/2). If/3-|-7r/2 < 7r(l/a-|-p), we can repeat the above 
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procedure: choose x in the sector arg(a;) G (/9,/9 + 7 r/ 2 ) and rotate the contour of integration by an angle 
[3 clockwise. Eventually we will cover the whole sector arg(a;) G (— 7 r/ 2 , l/a + p). Finally, we extend 
/^(x) into the sector arg(a;) G (—1/a — p,'n'/2) by the symmetry principle (the conjugate of fj^{x) is 
fj^{x)). This ends the proof of part (i). 

Let us now prove part (ii). Everywhere in this proof we will denote by Ai some positive constants 
which can depend only on (a, p). We choose any {3 G (0,7r(l/a — p — 1/8)) (if this interval is empty we 
take = 0). We set x = in (42) and rewrite that equation in the form 

POO 

/ e-Xe-'^z)dz. (43) 

Jo 

Note that (3 < 7 r(l/a — p), so the second term in (42) vanishes. Formula (43) is valid in the half-plane 
Re{y) > 0. The function p{u) is analytic in the sector | arg(M)| < 7 r(l/a — p), which includes the sector 
I arg(M)| < 7 r(l/a — p — 1 / 8 ). Given the behavior of piu) as u —)■ O'*" or u —)■ -|-oo given in (38), we see 
that there must exist a constant Ai such that 

|/i(n)| < Ai min(|M|", |m|“"^), (44) 

for all u in the sector | arg(M)| < 7r(l/a —p— 1/8). Note that Re(j/) > \y\/2 in the sector | arg(?/)| < 7r/4, 
thus for all y in the sector | arg(p)| < 7r/4 we have 

POO 

|/^(e"‘^l/)| < kli / e-"^l^l/2j^in(z“,W“^)dz. (45) 

Jo 

We leave it as an exercise to check that the above estimate implies that there exists a constant A 2 such 
that 

for I arg(p)| < 7 r/ 4 . Since the sector | arg(a;)| < 7 r(l/a — p-|- 1/8) can be covered by hnitely many sectors 
of angle 7r/2, we see have proved that there exist some constant ^3 such that for all x in the sector 
I arg(x)| < 7 r(l/Q; — p -|- 1 / 8 ) we have 

\f^{x)\ < A^mm{\xr-\\x\-^-^). 

Next, we take any [3 G (7r(l/a — p + 1/8), 7r(l/a — p -|- 1)) and repeat the same procedure. We set 
X = A^y in (42) and rewrite that equation in the form 

POO 

fx{A^y) = z)dz — 27ri x Res(p(M) : u = (46) 

Jo 

The integral term is estimated as above, and the exponential term is estimated as 

\e-(u-y^y\ = exp(-Re(e-"‘(^/“-^)+^^+‘”*^(*'))|p|) = exp(-cos(/? - 7r(l/a - p))|p|). 

This ends the proof of item (ii). 

The ideas for the above proof of item (ii) are illustrated in Figure 1. The contour of integration 
in (42) is over the red line With this contour of integration, the integral representation (42) 

is valid in the half-plane AOE, and we obtain asymptotic estimates of fx{x) in the sector BOD. The 
asymptotic estimates can be made uniform, as long as the contour of integration does not pass 

too close to the pole of p{u) as u = U-. 
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Figure 1: Illustration to the proof of Lemma 6. 


Let us now prove part (iii). Setting /3 = 1/a in the equation (46) and using formula (37) we obtain 


e"/“/x(e"/“a;) = -i 


Next, using (29) we compute 


.>S* 2 (ap)_ ^„nip I I o_; „/oN i f za; / —7ri/o 


exp(-a:e"^ + 7riapp/2 + 37rip/2) + / iJi{e-^^'^z)dz, (47) 

Jo 


n{e = — sin{7rap)z^^^^e ’^‘^/^52(a + ap/2 — ialn(2;)/(27r))52(l + a + ap/2 + ialn(2;)/(27r)) 

71 

The functional equation (9) gives us the identity 

S' 2 (a + ap/2 — ia ln(z)/(27r)) = 2sin(7r(l + p/2 — iln(^)/(27r)))52(l + a + ap/2 — ia ln(2;)/(27r)) 

= i {f'^Pl‘^z^l‘^ - ^2(1 + a + ap/2 - ia ln(^)/(27r)) 

and simplifying the result we obtain 


^(g-W«^) = -sin(7rap);2“^/2-i/2^g^ip^ _ ^ ^ ^ + ia ln(; 2 )/( 27 ri))| 2 . 

TT 

Combining the above result with (47) and (11) we get 

exp(—+ 7riapp/2 + 37rip/2) — — sin(7rap)(G(a;) + e’^‘^G'(a;)) 
va 71 


a 


= ^^^ 2(1 + ap) X — isin(7rp) exp(—+ 7riapp/2 + 37rip/2) 


(G(a;) + e"^G'(a;)) 


a sin(7rap) 
27rS'2(l + ap) 


Using formulas (9) and (A.7) we check that 


asin( 7 rap) \/a ^ , . 

= —h2(-ap), 


27rS'2(l + ap) dvr 
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and we leave it as an exercise to verify that 

—isin(7rp) exp(—+ TTiapp/2 + 37rip/2) = K{x) + 


where 

k{x) := sin(xsin(7rp) + 7rp(l — ap)/2). 

The above fonr formulas combined with the dehnition of F{x) (see formula (12)) imply the desired iden¬ 
tity (41). □ 


Remark 2. The proof of part (i) of Lemma 6 shows that fj^{x) can in fact be extended to an analytic 
function on the whole of IZ, though this result will not be used anywhere in this paper. Also, note that 
(41) implies the following results, 

^52(1 + ap)F{x) = f exp{e-^^P{y - x))/x(e"/“|/)d 2 / 

V ® Jo 

sm(7rp) ^ ' 

which establish a direct link between the spectral theory of stable processes on the half-line and the 
Wiener-Hopf factorization theory. 

Let us consider what happens to the asymptotic behavior of the function f^{x) (as |x| —)■ -|-oo) as 
we increase |arg(x)|. Lemma 6(ii) tells us that there is a transition from power-type decay to 

exponential growth exp(—^|g|) which occurs at the critical level | arg(x)| = a := 7r(l/a — p + 1/2). For 
values of | arg(x)| < a we have ^ > 0, thus we have a power-type decay; precisely at the critical level 
|arg(x)| = a we have ^ = 0 which results in a bounded oscillatory behavior; hnally, when |arg(x)| 
exceeds the critical level cr we have ^ G (—1,0) and fx{x) exponentially increases and oscillates. It is 
easy to check that for all (a,p) G ^ we have a > 7i/{2a). Using this result we rewrite the upper bound 
in Lemma 6(ii) in the following less informative (but more useful) form: 

Corollary 2. There exist constants C = C{a, p) > 0 and 9 = 9{a, p) > 0 such that for all q in the sector 
I arg(g)| < 7i/2 + 9 we have 

l/x(F‘*)l (48) 

and for all q in the sector 12 + 9 < \ arg(g)| < n we have 

{C min(|g|^“^/", if p < 1/2, 

|/x(?'/“)| < C{\q\P-^/- + 1), if p = 1/2, (49) 

\c{\q\P-^/^ + exp(|g|i/“)) if p > 1/2. 


The corresponding results for /x(x) can be obtained by changing p p. 

The next corollary will play a crucial role in the proof of Theorem l(ii) and Theorem 2(ii). 
Corollary 3. Assume that {a,p) G A. Then for all x,y > 0 we have 


r‘min(x,f/) o 

Im [e2"/“/x((x - z)e^'^/^)fx{{y - z)e-‘/“)] dz = -F{x)F{y). 


a 


(50) 
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Proof. First of all, we use formulas (9) and (A.7) to check that 


52(1 + ap)S2{l + ap) = 


2 sin(7rp) 


Next, we use identity (41) and compute 

Im [e"/“/x((x - z)e"/")/y((!/ - z)e’'l‘-)] 

^ ' \(F{x-z)+ €''•<'F'{x-z)){F(y-v)+e'''''F'(y - z)) 


■ ( 

asm(7rp) 

2 r , - 

— F{x — z)F {y — z) + F {x — z)F{y — z) 
a L 


2 d 
a d^; - 


F[x — z)F{y — x) 


To hnish the proof one needs to integrate the above expression in z G (0, min(a;, y)) and use the fact that 
F(0+) = -F(0+) = 0 (see Lemma 4). □ 


Proof of Theorem 1, part (i): Let us denote by mt{y) the probability density of the random variable 
—X_t when Xq = 0 (the existence of this density is known - see [4, 5, 14], though it also follows easily 
from the absolute convergence of the integral in formula (51) below). As we have discussed in Section 
2.1, the scaling property of stable processes implies that the probability density function of 
given by Conditioning on the exponential random variable e(g) we arrive at the Laplace 

transform identity 

e~‘^^mt{y)dt = 

We write down mt{y) as the inverse Laplace transform and then pass to Hankel’s contour (via Lemma 2 
and Corollary 2; here we need a>lorp>l/2): 

1 /■ 1 

My) = ^ dg. (si) 

Jm ^ Jo 

Next we compute 



Px(,To >t) = Po(-Xj <x)= mt{y)dy 

Jo 

Im [e"/“/x(|/g^/"e"/“)] dgj dg, 

where the application of Fubini’s Theorem is justihed due to Corollary 2. The integral in curly brackets 
can be evaluated using formula (41): 

r Im [e"/“/x(|/g'/“e"/“)] dy = ^^ 2(1 + ap) sin(7rp) T F'M/^)dy = ^^^g-'/“F(xgi/“), 

Jo ^ y/a Jo Va 

where we have used the fact that F(0+) = 0 (Lemma 4) and identity (9). Formula (13) follows by 
combining the above two identities and changing the variable of integration g = A". □ 
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Proof of Theorem 1, part (ii): We recall our notation from Section 2.1: the function ht{x,y,z) 
denotes the joint probability density of random variables {Xt,2Lt) ^ stable process X started from 
Xq = X. We start with equation (24), which we reproduce here for convenience: 

h,{x,y,z) = ^ f (52) 

JiR 

According to Corollary 2, for hxed 0 < z < min(x, y) the integrand is bounded from above by C min(l, \q\~^) 
when q G iM, thus the integral in (52) converges absolutely. 

Now we plan to pass to Hankel’s contour of integration in (52). Corollary 2 tells us that in some 
sector 7r/2 < | arg(z)| < 7r/2 + 6 the integrand is bounded from above by C min(l, l*?!”^). In the sector 
7r/2 + 6^ < I arg(^)| < tt the term is exponentially decaying. When a > 1, the integrand is bounded 
by which decays to zero exponentially. When a < 1 and p = 1/2 the integrand is bounded 

by C'e'^*, which again decays to zero exponentially. Thus all conditions of Lemma 2 are satished and we 
obtain 

1 r°° 

ht{x, y, z) :=- Im - z)e-'^/^q^/^)fx{{x - 

^ Jo 

Assume hrst that x ^ y. We use formula (22) and apply Fubini’s Theorem to obtain 

l‘min(x,y) 

Pt{x,y)= ht{x,y,z)dz (53) 

Jo 

1 poo I pmm{x,y) | 

= Im [e2-/“/5^((|/ - z)e"'/“g'/“)/x((a: - z)e"/“gi/“)] dz \ 

To justify the application of Fubini’s Theorem, we use Corollary 2 and estimate the integrand as follows: 

|/x((9 - (54) 

< Cq^l“-'^e-*((y - + 2l((/))((x - + B(q)], 

where A{q) = S(g) = 1 if p = 1/2 and a < 1 and A{q) = B{q) = exp(max(x, ?/)|g|^/“) if a > 1. 
Considering these two cases separately (a < 1 and a > 1) one can check that the function in the 
right-hand side of (54) is integrable in (z, q) E {0 < z < min(a;,|/), g > 0} C as long as x ^ y. 

To hnish the proof, we compute the integral in curly brackets in (53) using Corollary 3 and change 
the variable of integration g = A“. To remove the restriction x ^ y we note that the right-hand side in 
(14) dehnes a continuous function of t,x,y > 0, and we dehne pt{x,x) by continuity. □ 


2.5 Proof of Theorem 2 

Recall that we have dehned ( := |min(l, 1/a). Let us dehne Sr to be the “shift operator”: Sru{x) = 
u{r + x), and let Cu{z) := u(x)e~^^dx denote the Laplace transform of u. Note that u G "Ha implies 
SrU G "Ha for all r > 0. 

Lemma 7. Let u G Jia- Then CSru{z) is an entire function and for every small e > 0 there exists C = 
C{u,a, e) such that for allr > 0 we have \CSru{z)\ < C min(l, |^|“^) in the sector \ arg( 2 ;)| < ^-|-7r/2 —e. 
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Proof. The fact that CSru{z) is an entire fnnction follows from Dehnition 1, which implies that u{x) is 
bonnded for small x = 0 and decays faster than exp(—(5xln(x)) for large x. 

Let ns £x /? G (—C;C)- By the same estimate as above, we can rotate the contonr of integration and 
change the variable of integration x = to obtain 

poo 

CSru{z) = + P^y)dy. (55) 

Jo 

According to Dehnition 1, there exists 6{u,a,e) > 0 and C{u,a,e) > 0 snch that |M(r + e^^y)\ < 
C'min(l, for all |/ > 0. Using the estimate Re(e^^z) > sin(e)|z| in the sector | arg(e'^z)| < 7r/2 — e 
we obtain (for all z in this sector) 


\CSru{z)\<C < 


C 


sm e R 


Fnrthermore, 


CSrU{z)\ < C 



y ^y<ly. 


Since {3 can be chosen arbitrary in the interval (—C, 0) we obtain the desired resnlt. 


□ 


Next, for u G PLa we dehne 

poo 

Eu{\) = / fx{Xx)u{x)dx. (56) 

Jo 

Lemma 8. Let {a,p) G A and u G PLa- 

(i) There exists e > 0 such that for every r > 0 the function X i—)■ ESru{X) is analytic in the sector 
I arg(A)| < 7r/a + e where it satisfies 

|Hi5j,m(A)| < C X min(|A|“^“^, |A|“^), (57) 

for some constant C = C{u, a, p,e) > 0. 

(a) The function X G (0, oo) i—>• nn(A) satisfies 

|nM(A)| < C X min(l,A“^), (58) 

for some constant C = C{u, a, p) > 0. 

Proof. First we will prove part (i) in the special case r = 0 (so that S^-u = u). Recall that the Wiener- 
Hopf factor 0(z) (given by (26)) is the Laplace transform of fxi^) and that (f^z) can be extended to a 
meromorphic fnnction on 77 (the Riemann snrface of the logarithm fnnction). Moreover, according to 
(26) and properties (i) and (ii) of the donble sine fnnction on page 4, we see that the fnnction 0(x) is 
analytic in the sector | arg(z)| < 7r(l/a + p) and satishes \(t){z)\ < C'min(l, \z\~'^p) in any smaller sector 
I arg(x)| < 7r(l/a + p) — e for some C = C{a, p, e). 

Let ns take u G PLa- Using Plancherel’s Theorem for Laplace transform we rewrite (56) as follows 

Eu{X) = [ (j){z/X)Cu{-z)dz. 

27riA 
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If a < 1 we choose any e G ( 0 , 7 r/)), and if a > 1 we choose any e G (0, |(1 — l/a)). Let us dehne 
[3 = 7 r /2 — ^ + e (where C = f niin(l, l/a)). Applying Lemma 1, the upper bound for cj) and the 
properties of Cu{—z) described in Lemma 7 we can deform the contour of integration in the above 
equation iM i—>• 7 , where 7 = {e~^^oo, 0) U (0, e^^oo), and obtain 


im(A) = 


27riA 


(j){z/X)Cu{—z)dz. 


Rewriting this expression as two integrals and changing the variable of integration we arrive at 


Sn(A)=/i(A) + / 2 (A) 

i/3 poo -i/3 poo 

= —— / (j){we^^/X)Cu{—we^^)dw -— / / X)Cu{—we~^^)dw. 

2vriA Jo 27riA Jq 


(59) 


Since is analytic in the sector | arg( 2 :)| < 7r(l/Q; + p), the functions /i(A) and / 2 (A) are analytic in 
the sector | arg(A)| < nil/a + p) — /3. Considering separately the two cases a < 1 and a > 1 we check 
that ^ := 7rp — /d > 0, thus the function Eu{X) is analytic in the sector | arg(A)| < vr/a + ^. 

According to Lemma 7, we can estimate |£m(—< C min(l, |w|“^). This result and the bound 
< C min(l, \w\~'^^) combined with Lemma 3 give us the upper bound (57). 

The proof of the case r > 0 follows exactly the same steps, since the function Vr{x) = Sru{x) also 
belongs to Tia- Lemma 7 then tells us that the upper bound in (57) is uniform in r. 

The proof of part (ii) follows very similar steps. We choose e and (3 in the same way as above. Let 
f{z) := \/2/Tr F{x)e~^^dx. As we have established in formula (19), the function f{z) can be analyt¬ 

ically continued to a meromorphic function on 77. Moreover, f{z) is analytic in the sector | arg(z)| < up 
and in any smaller sector |arg(z)| < vrp — 62 it satishes an upper bound \ f{z)\ < Cmin(l, \z\~^~‘^p), 
which follows from (10) and (19). 

Note that e~^F{x) and e^^u{x) are in //^(M"*"). We express IIm using Plancherel’s theorem for Laplace 
transform as follows 


Uu{X) = / {e-^^F{Xx)){e^^u{x))dx = 


2mX 


f{f)Cu{-z)dz. 


'A+iR 


The contour of integration is deformed in the same way as in (59) and we obtain 


r°° .r, .r. P~^P 

nM(A) = / f{wF^/X)Cu{-wF^)dw - 

V 0 

Applying Lemma 3 and the above mentioned estimate \f{we~' 
result (58). 


f{we /X)Cu{—we ^^)dw. 

^)\ < C min(l, we get the desired 

□ 



Lemma 9. Let a and p he as in Theorem 2. Then 11 : L^(]R’'') 1— )■ //^(M"*") is a bounded self-adjoint 
operator and II : TLa 1 —)■ L^(]R+) is a symmetric operator such that TlTLa is dense in L^(M+). 

Proof. The fact that both operators 11 and II are symmetric follows easily from their dehnition (15). Let 
us denote by M. the Mellin transform operator: 


1 

Muiz) = . _ / u(x)x~^^‘^^^^dx. 


(60) 
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It is well-known that the Mellin transform operator is an isometry between and Since 

the operator 11 is a Mellin convolntion, Theorem 3(ii) and a standard compntation leads to identity 

M.Tlu{z) = MF{z) X Mu{—z). 

Formnla (20) and location of poles/zeros of the donble sine fnnction imply that M.F{z) has no poles on 
the real line, and the asymptotic relation (10) tells ns that \AiF{z) \ exp((p — l/2)7r|2;|) —)■ 1 as 2 ; —)■ oo. 
Thns for p >1/2 the fnnction AiF{z) is bonnded and it indnces a bonnded mnltiplication operator on 
L^(M), thns the operator Hu = M.~^[M.F x M~^u] is a bonnded (self-adjoint) operator on L^(M+). 

The fact that fl maps Ha into L^(M+) follows from Lemma 8(ii). Let ns prove that the image 
tVHa is dense in Assnme to the contrary that IlHa is not dense in Then there mnst 

exist a nonzero fnnction w G L^(M+) such that tc(a:)IlM(a;)da: = 0 for all u G Ha- Consider any 
nonzero function u G Ha and dehne Ua{x) = u{ax) for a > 0. As we have mentioned in the paragraph 
following Dehnition 1, we have Ua G Ha for all a > 0. Let n(A) = nn(A). Then it is clear from the 
dehnition of fl that IlMa(A) = a“^n(a“^A), so that all functions v{a~^x) belong to TVHa and we have 
w{x)v{a~^x)dx = 0 for all a > 0. Applying Plancherel’s Theorem for Mellin transform we see that 

0 = / w{x)v{a~^x)dx = [ M.w{—z)M.v{z)a^^‘^^'^^dz, for all a > 0. 

Jo Jr 

Thus the Fourier transform of a function Aiw{—z)Aiv{z) G L^(M) is identically zero, therefore 

Mw(—z)Mv(z) = 0, for all z G M. (61) 

Next, due to the fact that v = Iln and the estimate (58) we see that M.v{z) is analytic in a strip 
|Im( 2 ;)| < 1/2, thus A4v(z) is nonzero almost everywhere in this strip. Equation (61) then implies that 
Mw(—z) = 0 for almost all ^ G M, therefore w = 0 and we arrive at a contradiction. □ 


Proof of Theorem 2: Given the result in Lemma 9, it remains to establish formulas (16), (17) and 
(18). Our plan is to prove (17) using similar method as in the proof of Theorem l(ii) and then to derive 
the remaining statements (16) and (18) as corollaries of (17). 

When p = 1/2 formula (17) follows at once from (14) and Fubini’s Theorem (whose application is 
justihed since the functions F and F are bounded in this case). Thus we will only consider the case 
when p > 1/2. 

We recall that the operator S is dehned by (56) and Sr is the “shift” operator from Lemma 7. Let 
us hx M G Ha- We begin by rewriting formula (22): 

px 

Ptu{x) = / kt{x,z)dz, (62) 

Jo 

where we have dehned 

POO 

kt{x,z):= ht{x,y,z)u{y)dy. 

J Z 

The Laplace transform of f h->■ kt{x^ z) can be computed by applying Fubini’s Theorem to the above 
formula and using (21), (23) and (56): 



e ‘^*kt{x,z)dt 


/ OO 

q~^Hq{x, y, z)u{y)dy = g^/""Vx((a: - 2;)g^/“)S5^M(g^/“) 


Kq{x,z). (63) 
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Let us denote w = {x — z)\q\^/°'. According to Corollary 2 and Lemma 8, the function q i—)■ Kg{x, z) 
is analytic in the sector | arg(g)| < tt + e and it is bounded in the sector | arg(^)| < vr by 

\Kg{x,z)\ < C{u,a,p) X min(t(;"^“^, x min(|g|^“^'^“, (64) 

For fixed z and x the upper bound (64) implies the estimate: 

|A'g(a;, x)| < C(M,Q:,p,a;, x) X min(l, |g|“^), |arg(g)|<7r. (65) 

At the same time, using the trivial bounds min(t(;“^“^, and mindgl^"^/", \q\~^^^) < 

we obtain from (64) the following result: 

|A’g(a;, x)| < C(m, Q:,p) X |a; — |arg(g)| = 7r. (66) 

Next, we fix x,z satisfying 0 < x < x, use the upper bound (65) and express kt{x,z) as the inverse 
Laplace transform 

kt{x,z) = — Kg{x,z)e'^^dq. 

Jm 

Next, we use the same upper bound (65) and Lemma 2 and pass to Hankel’s contour of integration in 
the above integral: 

I f°° 

kt{x,z) = -/ Im[iL„g,ri(x, x)]e“^Mg. 

^ Jo 

Combining this result with formula (62) and applying Fubini’s Theorem we obtain 


1 r°° 

Ptu{x) = — / 

71" .In 


Im[iFge’ri(x, z)]dz 




(67) 


The use of Fubini’s Theorem in the previous step is justihed since the function (x, g) iCgg,ri(x, x)e“'^* 
is absolutely integrable on (0,x) x (0,oo) - this follows easily from the upper bound (66). 

Next, we rewrite the integral in square brackets in (67) as follows 


Im[iCqg,ri(x, z)]dz 


( 68 ) 


= -g2/“-i 


= -g2/“-i 


X /*oo 



Im - z)q^'V''‘)fx((y “ dt/dz 

'0 Jz 

r*oo pmm{x,y) 

Im <izu{y)dy 


'0 ^0 

poo 


= -gi/"-^ / F(xg^/")F(j/g^/“)M(j/)d|/ =-g^/“-iF(xg^/“)nM(g^/"). 


a 


a 


In the second step we have again applied Fubini’s Theorem: to justify its use note that according to 
Corollary 2 we have 


\fx{(x - z)q^'V‘'“)fx((y - z)q'/"e"/’^)u(y)\ < C X [x - X \y - z^z-^e^ X |u( 9 )| 


for some constant C = C{a,p,q), and the function in the right-hand side of the above equation is 
integrable over the region {{z,y) G : 0 < x < min(x,g)}. In the third step of (68) we have applied 
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Corollary 3. Combining formulas (67) and (68) and changing the variable of integration g = A“ we obtain 
the desired result (17). 

Formula (16) follows from (17) by taking the limit t —)■ O’*": Then Ptu{x) —)■ u{x) for all x > 0 (since 
u is continuous and a Levy process X killed on the first exit from (0, cx)) is continuous in probability), 
while the right-hand side of (17) converges to IHIm. 

Let us now prove formula (18). We denote by (u,u) the inner product in We start with 

v,w G Tia and use formulas (16) and (17) and the fact that operators If and If are symmetric to compute 

{Ptflw, v) = (nw, Ptv) = (Liu;, ne-^“*nu) = (nnu;, e-^“*nu) = (w, e-^“*nu) = (ne-*^“u;, v). (69) 

In the last step we have used the fact that G Tia- Since (69) is valid for all v G Pa and Pa 

is dense in L^(M’'') we conclude that Ptllw = Therefore, flPiLIrc = Lllle^^^'^w = as 

desired. □ 


Proof of Corollary 1: Part (ii) of Corollary 1 was established above (see the discussion following (69)). 
To prove part (i), we start with u G Pa, v G and denote flu = w. According to (16) we have 

flic = u. We use formula (18) and compute 

{PtUv, w) = {PtUv, flu) = (u, UPtflu) = (u, e-*^“u) = (e-*^“u, Hu;) = (ne-*^“u, w). (70) 

Since w ranges over 11770,, and the latter is a dense subset of L^(M+) due to Theorem 2(i), we see that 
(70) implies that for all v G //^(M"*") we have PtHv = Application of Fubini’s Theorem (note 

that ux{x) = F{Xx) is bounded when p > 1/2) gives 

poo poo poo 

Ptllv{x) = / / pt{x,y)F{Xy)v{X)dXdy = / PtUx{x)v{X)dX. 

Jo Jo Jo 

Since v G is arbitrary, we have PtUx{x) = e~^^°‘ux{x) for almost all x, A > 0. By continuity, this 

relation holds for all x, A > 0, as desired. □ 


3 Doney classes 

Doney [6] has introduced classes of stable processes for which the Wiener-Hopf factors can be computed 
explicitly in terms of q-Pochhammer symbols. Let fc, / G Z and let us define by Ck^i the class of stable 
processes with parameters (a, p) satisfying 


ap = I — ka. (71) 

These classes include the spectrally-positive (Co,i) and the spectrally-negative processes (C_i _i). Note 
that X G if and only if X G Let us define the q-Pochhammer symbol 


(a;g)n 


' n—1 

n(l-ag^'), n>0, 

i=o 

|n| 

n(l - ag ■^) ^ n< 0 , 
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and (a; g)o = 1. 

The main result of Doney [6] is the following formula for the Wiener-Hopf factor: for a process 
X G Ck,i we have 

((_l)«+l^ag7ria(l-fc).^)^ 




(72) 


where q := and q := e“^. We note that this formula follows from our general result (Theorem 
5) and formula (A.8). It is not surprising that all expressions involving the eigenfunctions F{x) also 
simplify considerably for Doney classes Ck,i- 

Proposition 2. Assume that X G Ckj- 

(i) The function G(x) (which is related to F{x) through (12)) is given by 




^■q) 

G{x)= I —- ry—J^dz. 


(73) 


(a) ForRe{z) > max(0, cos(7rp)) 


Ja ((-l)'z“e 

e-^^F{x)dx = ^S2{ap) ^ 


a „iTia{k+2). 




((-l)^Ze-"h«;g)_ 


/ + 1 


(74) 


(Hi) Assume that p > 1/2 and Re( 2 ;) G (—op, 0). If I > 0 then 

poo -I k-\-l I 

/ x^~^F{x)dx = -(—TT 2sin(7r(^ + (j — l)a)) TT 


while if I < 0 we have 


/ a;^“^F(a:)dx =-(—TT2sin(7r(2: + j — 1 )/q;) TT 


2sin(7r(z - j)/a)’ 


|fc+i| 


2sin(7r(2; — jo)) ’ 


(75) 


(76) 


Proof. The proof follows immediately from (11), (19), (20) and formulas (A.8) and (A.9). 


□ 


3.1 Spectrally one-sided processes 

The spectrally-negative (spectrally-positive) processes belong to classes Co,i (respectively, C_i _i). All 
the formulas given in Proposition 2 for general Doney classes Ck^i remain valid, except that in the 
spectrally-negative case op = 1 we should take G{x) = 0. Let us consider this case in more detail. 

In the spectrally-negative case it is well-known that the random variable We(i) has exponential dis¬ 
tribution with parameter one, so that 


fxi^) = e a; > 0. 

Following the same analysis as before (or simply using formula (41)), we would obtain the co-eigenfunction 

F{x) = sin(xsin(7r/a)). 
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The eigenfunction is given by (12) and (73): 


F{x) = e*“s(7r/a) sin(xsin(7r/Q;) + 7r(2 — q;)/(2q;)) 


a 


sin ttq; 


/o 1 + 2 cos(7rQ;)M" + 


4 Conclusion 

Our methods for establishing Theorems 1, 2 and 3 are fundamentally based on the Wiener-Hopf fac¬ 
torization (Theorem 4), which comes from the fluctuation theory of Levy processes. We hope that the 
methods developed in this paper will help to study spectral properties of semigroups of more general 
Levy processes on a half-line (see [16] for some results in this direction). Of course, the case of stable 
processes is a very special one, as stable processes lie in the intersection of the class of Levy processes and 
the class of self-similar Markov processes. We would like to mention that the spectral properties of the 
semigroup of general positive self-similar Markov processes (which include stable processes killed on the 
first exit from (0, oo)) are studied in [21] using a different approach (based on Lamperti transformation 
and Mellin transform techniques). 

Finally, we would like to emphasize the importance of the double sine function S 2 {z] a) in the study 
of stable processes. This function appears in the definition (12) of the eigenfunctions, in formulas (19) 
and (20) that give the Laplace and Mellin transform of the eigenfunctions and in formula (26), which 
gives the Wiener-Hopf factors of stable processes. So far, apart from some number-theoretic applications 
[12, 26], this curious special function has appeared mostly in the Physics literature [8, 22, 24, 29], where 
it is used in studying quantum topology and cluster algebras. It is an interesting question whether this 
appearance of the double sine function in our study of stable processes is simply a coincidence or there 
is indeed a deeper connection between stable processes on the half-line and quantum topology and/or 
cluster algebras. 
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Appendix A The double sine function 


In a series of papers Koyama and Knrokawa (see [11, 12], for example) have stndied the mnltiple sine 
fnnctions Sr{z; w) where ^ G C, r G N and w = (tCi,..., Wr), w* G C \ (— cxd, 0]. These are homogeneons 
fnnctions (in the sense that Sr{cz-, cw) = Sr{z; w) for any c > 0) that satisfy certain fnnctional eqnations. 
Considering the case r = 2, the fnnetion S 2 {z;w) satisfies 


S'2(2; + tCj;w) 


>5'2(^; w) 

2 sm{7iz/wi) ’ 


*g{1,2}. 


It can be shown [28] that 5*2 ( 2 ;; w) is a nniqne meromorphic fnnetion of z which satishes the above 
fnnctional eqnations and has valne one a.t z = {wi + W2)/2. Dne to the homogeneity of S 2 {z; ( 101 , 102 )) 
we can set Wi = 1, and from now on we will write simply S 2 {z;a) = S 2 {z; (1,0;)). Note that the same 
homogeneity property implies an identity 


S 2 {z;a) = S 2 {z/a;l/a). (A.l) 

As the donble sine fnnetion is becoming better known and more widely nsed in Mathematics and 
Physics literatnre, it seems that most anthors introdnee a new notation for this object, which makes it 
very hard to navigate the literatnre and to find relevant resnlts. Below we provide a short snmmary of 
varions notations nsed for the donble sine fnnetion S 2 {z;a): 

(i) The donble sine fnnetion Sb{z) of Ponsot and Teschner [22]: 

S2{z;a) = S^{z/^)~\ 


(ii) The hyperbolic gamma fnnetion Ghyp(wi,tC 2 ; z) of Rnijsenaars [24, 25]: 

S 2 {z-, a) = Ghyp(l, a; i(l/2 + a/2 - 


(iii) The hyperbolic gamma fnnetion Vh(z-,wi,W 2 ) of van de Bnlt [28]: 

S2{z-,a) = Th{z-, l,a)"^ 
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Figure 2: The roots (blue crosses) and the poles (red circles) of the double sine function 5*2 (^; a) 


(iv) The double gamma function G{z] a) of Barnes [1]: 


S2{z-,a) = 

^ ^ ^ ’ G{l + a- z]a) 


(A.2) 


(iv) The 7 -function (see [29] and the references therein): assume that Im(a) > 0 and dehne q := e'^™, 
q := and 

where (a,g)oo := nA:>o(^ “ q-Pochhammer symbol. The y-function is related to the 

double sine function through the identity 


S2{z;a) = g-’^i(2^-l-a)^/(8a)-|-7ri(o-|-l/a)/24 ^ 


7(z;a)' 


(A.4) 


(iv) The quantum dilogarithm function <F^(z) (see [8] and the references therein) is related to the double 
sine function S 2 {z] a) through the above formula (A.4) and the identity 


■-y{z; a) = <F"(7ri(l + a — 2z)) 


-1 


For hxed a G C \ (—cxo,0], the double sine function is meromorphic in ^-variable and it admits the 
product representation of the form 


S 2 {z-a) = e^G)+B(o.)z+C(a)z^ 


P{—z/{m + an)) 

z — l — a^^^^ P{z/{m + 1 + a{n + 1))) 

m>0n>0 ^ ^ 


nn 


(A.6) 


where P{z) := (1 — z)exp{z + z‘^/2) and the prime in the second product means that the term corre¬ 
sponding to m = n = 0 is omitted. The explicit expressions for A{a), B{a) and G{a) can be found in 
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[26, Theorem 6]. It is clear from (A.5) that S 2 {z]a) has zeros at points {—m — an : m > 0,n > 0} 
and poles at points {m + an : m > l,n > 1}. Moreover, the pole at z = a + 1 is simple (other poles 
may have mnltiplicity greater than one when a is rational). Note that if Im(a) > 0 and if we look at 
the lattice {m + an, m,n E Z}, then the poles (zeros) of S 2 {z;a) he in the hrst (respectively, third) 
qnadrant of this lattice (see Figure 2). 

Besides the value S'2(l/2 + a/2;a) = 1, the following values are known explicitly [12]: 

S' 2 (l; a) = 1/S'2 (q;; a) = \/a, S' 2 (l/ 2 ; a) = S' 2 (q;/ 2 ; a) = a/2. (A.6) 

The above result (A.2) implies a very useful reflection formula 

S 2 {z-, a)S 2 {l + a — z] a) = 1. (A.7) 

Let us denote q = and q = e~^. Using (A.3), (A.4), (A.7) and the identity (a; g)oo = 
{aq^] q)oo{(^] <i)n we can derive a useful result 


S' 2 (l /2 + a/2 + (m — na)/2 + a)S' 2 (l /2 + a/2 + (m — na)/2 — z] a) 

a i \m—l^— 27 r\z-\- 7 ria(l—n). 
((_l)n-lg- 27 rn/a+ 7 ri(m-l)/a. ' 


(A.8) 


= e 


When m,n E Z+, we can apply repeatedly the functional equations (9) and prove the following result 


S2{.z]a) 

S 2 {z + m — na] a) 


' l^mn TT 2 sin(7r(^ + j — l)/a) TT - ^ 

^A^A2sin(7r(2-ja)) 


(A.9) 


A similar expression when —m, —n E Z"*" can be easily obtained from the above formula by changing 
variables z + m — na = w (or by using transformation (A.l) or the reflection formula (A.7)). 

An asymptotic behavior of the double sine function as the imaginary part of the argument increases 
is described by the following equation 

|S' 2 (c + iy; a)| = + o(l)), y —)■ oo. (A.10) 


See formula 135 in [22]. It is easy to prove that this result holds uniformly for c on compact subsets of 
M. We note that this result implies (10). 

Finally, the “b-beta integral” of Ponsot-Teschner (Lemma 15 in [22], see also formula 7 in [29]) implies 
the following result: For a > 0, 0 < 5 < (1 + a)/2 and s E (—6, b) we have 


x* ^|S'2(l/2 + a/2 + 6 + ialn(a;)/(27r);a)Fda; 


27r 52(25; a) 

y/a ^ 2(5 + s; a)S 2 {b — s; a) 


(A.ll) 
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